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We report extensive simulations of liquid supercooled 
states for a simple three-sites molecular model, introduced by 
Lewis and Wahnstrom [L. J. Lewis and G. Wahnstrom, Phys. 
Rev. E 50, 3865 (1994)] to mimic the behavior of ortho- 
terphenyl. The large system size and the long simulation 
length allow to calculate very precisely — in a large q— vector 
range — self and collective correlation functions, providing 
a clean and simple reference model for theoretical descrip- 
tions of molecular liquids in supercooled states. The time 
and wavevector dependence of the site-site correlation func- 
tions are compared with detailed predictions based on ideal 
mode-coupling theory, neglecting the molecular constraints. 
Except for the wavevector region where the dynamics is con- 
trolled by the center of mass (around 9 nm" 1 ), the theoretical 
predictions compare very well with the simulation data. 



I. INTRODUCTION 

In the last decade, increased computational resources 
have been used to tackle the study of the onset of glassy 
dynamics [ffl-B, one of the most interesting open prob- 
lems in the physics of liquids. Supercooled liquids are 
indeed characterized by an extreme T-dependence of the 
structural times, which covers more than 15 orders of 
magnitude in a small temperature interval. Present com- 
puter facilities allow to follow the change in the structure 
and in the dynamics of the system (in equilibrium) over a 
T range where the diffusivity changes more than 5 order 
of magnitudes. At the same time, computer simulations 
are starting to provide a detailed picture of the structure 
of the potential energy landscape probed in supercooled 
states l^-H. The calculated trajectories offer an out- 
standing possibility for studying the origin of the slowing 
down of the dynamics, and call for a careful comparison 
between theories and "exact" numerical calculations. 

"Exact" results calculated for model-system and theo- 
retical predictions based on mode coupling theory (MCT) 
[ plpl have been reported for several atomic liquids, - 
encompassing hard spheres, soft spheres, Lennard-Jones 
PL and recently silica ]l7| |. In the case of molecular liq- 
uids, where the possibility of comparing with experiments 
makes the effort even more valuable, molecular dynam- 
ics (MD) simulations extending to the 100 ns timescale 
are now possible. Detailed comparisons between theories 
and simulations appears as well |p^|-^l||. 

MD simulations have often been used to generate ref- 
erence systems for validating theories of liquids, starting 



from the pioneering work of Alder and Wainwright |22| ] 
on hard spheres. A very interesting model for a molec- 
ular liquid was designed by Lewis and Wahnstrom (LW) 
p3[ by gluing in a rigid molecule three identical Lennard 
Jones (LJ) atoms. The shape of the molecule (an isoscele 
triangle) and the LJ parameters were chosen to mimic as 
close as possible one of the most studied glass-forming, 
liquid |J-|l| ort/io-tcrphenyl (OTP). This model pro- 
vides a clean bridge between studies based on atomic LJ 
potentials and more complex molecular potentials. From 
a theoretical point of view, such a model is also very 
valuable, since the LW model lends itself to several ap- 
proximation of increasing complexity. For example, one 
could consider the LW system as a liquid of identical 
atoms, correlated in space according to the site-site cor- 
relation function (without differentiating between intra 
and intermolecular sites), or one can improve by includ- 
ing the molecular correlation via the triplet site-site-site 
correlations. Solution of the exact molecular model, im- 
plementing a site-site description [ [32] ] which accounts for 
the intramolecular constraints or an expansion in gener- 
alized spherical harmonics |3^,|34| is also foreseeable. 

The LW potential was introduced in Ref. |2^], which 
reports a study of the slow dynamics of the LW liquid. 
The major properties of such models were discussed both 
for self and collective properties. In following papers [ p5[ , 
LW made also some preliminary attempts to describe the 
intermediate time dynamics using MCT predictions. Fa- 
vored by the increased computational power, we have 
decided to revisit such a model and — by extending 
the simulation time by a factor 20 and the number of 
molecules by a factor 30 — to calculate detailed proper- 
ties of the dynamics in supercooled states, like the non 
ergodicity parameters for both self, collective and orien- 
tational properties to be used as reference system to be 
compared with theoretical predictions. We present such 
data in this article. We also solve the MCT equations for 
such a model using as input the site-site structure factor 
and the site-site-site triplet correlation functions calcu- 
lated from the simulated trajectories. The comparison 
between the molecular dynamics data and the MCT pre- 
dictions is very valuable in assessing the role of the center 
of mass dynamics (which is neglected in the chosen MCT 
approach), the role of the triplet correlations ]|(||l7j and, 
in general, the ability of ideal MCT to capture the dy- 
namics of molecular liquids in weakly supercooled states. 
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II. THE MODEL AND NUMERICS 

The geometry of the LW molecule is a rigid three sites 
isoscele triangle hosting on each site a Lennard Jones 
atom. The site-site spherical potentials is 

l/(r) =4e((-) 12 -(-f) + A + Br (1) 
r r 

with e = 5.276 kJ/mol, a = 0.4828 nm, A = 0.4612 
kJ/mol and B = —0.3132 kJ/mol/nm. The parame- 
ters of the potential have been selected to reproduce bulk 
properties of the OTP molecule |2^] . The value of A and 
B has been selected to bring the potential and its first 
derivative to zero at r — 1.2616 nm. The resulting po- 
tential has a minimum at r = 0.5421 nm of depth —4.985 
kJ/mol. The length of the two short sides is 0.483 nm 
and of the long side 0.588 nm, forming an isoscele angle 
of 75 degrees. 

We have studied a system composed by N = 9261 
molecules for several state points, from T = 255 K to 
T = 346 K, as listed in Table |. The choice of a large 
system size was performed to avoid spurious oscillations 
in the correlation functions introduced by the periodic 
boundary conditions. The large size was also preferred 
to access information — here and in future analysis — 
on the small wavevector dynamical behavior. 

In this article we discuss in detail the auto correla- 
tion function for the collective and self density operators, 
both for center of mass (COM) and sites (s). The den- 
sity operators p (and their corresponding autocorrelation 
functions) are defined as j37j 

1 N * 

PcOM (q) = -7=* 22 e' q ' r C0M ; 

f com (q, *) = (pcom (q, t)p C oM (q, °)*> ( 2 ) 

, AT 3 

Ps (q) = — E e * ^ : F ° fo> *) = <P' (q, o) *) 

»=1 0=1 

(3) 

pf M (q)^ <q "' COM ; 
Foil (q. *) = (Pell (q. t)p s c l o{: (q. °)*> ( 4 ) 



We also present the /-dependent rotational correlators 
for the molecule symmetry axis u, defined as 

Ci(t) = (Pi[cos6(t)}) (6) 

where 6(t) = cos- l [{(i(t) ■ ^(0))/(/i(0) • /i(0)>] and P ; is 
the Z-order Legendre polynomial. The triplet site-site-site 
correlation function C3(k, pq) defined as 

(p s (k)p s (p)p s (q)) 
= 3NS a (k)Ss(p)S s (q)5q M+p (l + n 2 c 3 (k,pq)) (7) 

has been calculated evaluating the left hand side of Eq. [7]; 
n is the site number density. The evaluation of C3 has 
been performed at T = 266 K averaging over 1000 inde- 
pendent configurations, spanning a time interval of 75 ns. 
Where possible, up to 300 different triplets of wavevec- 
tors with the same moduli and relative angles have been 
averaged. We have calculated C3 for moduli of k,p and q 
less than 44.4 nm^ 1 with a mesh of 0.44 nm" 1 . 

III. MODE COUPLING APPROXIMATIONS 

In this article we present two theoretical calculations 
of the dynamics of the sites. The first approximation 
neglects completely the presence of the molecular con- 
straints and assumes that the site-site structure factor 
includes all requested information on the structure of 
the liquid. This strong approximation allows to apply 
in a straightforward manner the MCT equation for sim- 
ple liquids. According to MCT the time evolution of the 
normalized collective density correlation functions $> q (t) 
is given by 

$„(*) + i/,$ 9 (f) +n2$ B (t) + Jl* / dsm q (t-s)$ q (s) = 0. 

Jo 

(8) 

Here Q q = \fqv / ' S s {q), with v denoting the thermal ve- 
locity, is an effective phonon-dispersion law and v q = 
Viq 2 denotes a damping constant. The kernel m q is given 
as m q {t) = Fq({&k{t)}), where the mode-coupling func- 
tional T q is determined by the structure factor: 

mh}) = lf 70jtw fc /| q _ k |, (9) 



pff(q) = J*'**; Fr lf (q,t) = (pf f (<l,t)pff( q ,0y) 

(5) 

where r' COH and r* are the position of the COM and of 
the j-site of molecule?. Moreover, S COM (q) = F COM (q,0) 
and S s (q) = F s (q, 0) are the COM and site static struc- 
ture factors. 



Vq.k = S s (q)Ss(k)S s (\q - k|)4 [q • kc fc + q • (q - k) C , q _ k |] 2 . (10) 

where c q = (1 — S s (q)~ 1 )/n is the direct correlation func- 
tion in g-space. 

The second approximation is devised to retain some 
of the information of the molecular shape through the 
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inclusion of the triplet correlation function C3 |36|] , previ- 
ously assumed to be zero. The memory function is then 
given by 

71 

V^ = S s (q)S s (k)S s (\q-k\)^ 
[q-kc fe + q- (q-k)c| q _ k | - q 2 n c 3 (q, k, q - k)] 2 . (11) 

For simple atomic liquids, previous studies |3(| have 
shown that the approximation C3 = does not modify 
the MCT predictions, except for a small shift in the es- 
timate of the critical temperature. Recently, it has been 
shown [ p^t that in network forming liquids, the C3 con- 
tributions play a crucial role. Including the C3 contribu- 
tions could, in the present case, be relevant for describing 
the intramolecular site-site relations. For this reason, we 
compare the numerical data with both predictions. 

We numerically solved Eq. || on a grid of 300 equally 
spaced q values extending up to q = 120 rom -1 , imple- 
menting the efficient techniques described in Ref. [ p8| . 
The long time limit of & q (t) at the dynamical critical 
temperature — the so-called non ergodicity parameter - 
f q is obtained by an iterative solution of the bifurcation 
equation 

Y^r = mh})- (12) 

1 Jq 

Close to the critical point, 20000 iterations where re- 
quested to solve Eq.|l2| with a precision of 10 -15 . We 
have solved Eqs.|] for both approximations. The esti- 
mated critical temperatures are T = 152 K and T = 340 
K for the two approximations. Thus, neglecting com- 
pletely the molecular constraints (first approximation) 
strongly weakens the vertex function, resulting in a the- 
oretical critical temperature lower than the numerical 
T%id ■ Including the triplet correlation moves up the crit- 
ical temperature by a factor bigger than 2, resulting in 
an overestimation of the critical temperature, as always 
found in all previous MCT calculations. Notwithstand- 
ing the large difference in critical temperature the expo- 
nent parameter |l5| | is A = 0.70 in both approximations, 
corresponding to the exponents b — 0.64 and a = 0.32 

IV. MOLECULAR DYNAMICS RESULTS 
A. Static 

The site and COM static structure factors are shown 
in Fig. 0, as a function of temperature. The COM struc- 
ture factor does not show any appreciable change with T. 
Instead, the site structure factor increases significantly at 
the first peak position (corresponding in real space to the 
site-site nearest neighbor distance) , consistently with the 
behavior of simple liquids. The absence of any change in 



the COM structure factor shows that the slowing down of 
the dynamics in this model is controlled by the changes 
in orientational order. It is the ordering process in the 
orientational degrees of freedom which drives the slow- 
ing down of the dynamics. Any theoretical approach 
attempting to predict the slow-dynamics in this model 
must require orientational (or site-site) static informa- 
tion as input. The absence of any T-dependence in the 
COM structure factor, except for the very clear reduction 
of the compressibility on cooling (see below), suggest the 
possibility that T-independent peaks in the experimen- 
tally measured OTP structure factor are associated to 
COM features. Of course, the present model does not al- 
low a straightforward comparison with experimental data 
(as performed for example in Ref. |^|), since proton and 
carbon atoms are not explicitly taken into account. 

The large simulation box allows to precisely calculate 
the q — ► limit of the structure factor. As expected, 
we find that S COM (q -> 0) = S s (q -> 0)/3. The In- 
dependence of S s (q — > 0) is shown in the inset of Fig. 

& 

B. Self-Dynamics 

Fig. H shows the T dependence of the diffusion coef- 
ficient D, evaluated from the long time behavior of the 
mean square displacement. Compared to the LW original 
data, the range of T where equilibrium states has been 
simulated is larger and the precision in the determina- 
tion of the values of D is improved. The larger T-rangc 
shows that two regions, which differ in the T-dependence 
of D, are observed. Above T = 275 K, D(T) is con- 
sistent with a power law behavior in T — T^ D , with 
T C MD = (265 ± l)K and exponent 7 = 2.0 ± 0.1. Below 
T = 275K, D decreases on a slower pace, which could be 
approximated with an Arrhenius law. The two different 
behaviors reflect the dynamics differences associated to 
the region of validity of MCT (T > T^ D ) and to the 
region (T < T^. ID ) where hopping phenomena become 
relevant. The change in dynamical behavior shows up al- 
ready for \T — Tmd\/Tmd = 0-06, confirming the general 
expectation that hopping phenomena may mask the crit- 
ical dynamics close to the critical MCT temperature. As 
discussed at length in Ref. , the presence of hopping 
phenomena in this model makes impossible a clear-cut 
identification of the critical dynamics, without making 
use of information provided by the T— dependence of the 
a-relaxation phenomenon. Similar changes in dynami- 
cal behavior close to the critical MCT temperature have 
been recently reported in models of silica [ ftl| and wa- 
ter p2[| . It is worth noting that in the original work of 
LW, hopping phenomena in the orientational degrees of 
freedom were detected starting from T = 266 K. 

The temperature dependence of the site self-correlation 
function F^ el ^ is reported for a specific wavevector in 
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Fig. H Similar figures characterize the time dependence 
of F£ el f at different values of q. The change in dynam- 
ics above and below T% lD can be observed also in the 
so-called time-T superposition graph, where curves for 
the same correlator are shown as a function of a scaled 
time. Indeed, in the a-relaxation region the so-called 
time-temperature superposition principle states that it 
is possible to represent the T-dependence of an arbitrary 
correlator (p(t) on a single master curve, by rescaling the 
time via a single T-dependent timescale, i.e. 



0(i) = </>(t/r(T)) 



(13) 



The scaling time r(T) is conveniently defined as the 
time when the correlation function has decayed to an ar- 
bitrary chosen value, for example 1/e. Fig. || shows the 
time-T superposition for the COM self density correla- 
tion function both above and below (inset) Tmct- Only 
for T > 275 K, where D is well described by a power-law, 
the time-T superposition principle is well obeyed. Below 
275 K, deviations start to be noticeable and the plateau 
value starts to increase. 

To describe in a compact way the q dependence of 
the a-relaxation and to make contact with the theo- 
retical and experimental evaluations, we present in the 
next two figures the parameters which better describe 
the self density-density correlation functions according 
to the Von Schweidler law and the stretched exponential 
form. For a generic correlation function (frit), the Von 
Schweidler law, 

m =f e - h {1) (t/r) b + h (2) (t/r) 2b + 0{{t/ T f b ) (14) 

describes, including the next to leading order corrections, 
the departure from the plateau value f c in the early a- 
relaxation region. The amplitudes hru and /i( 2 ) strongly 
depend on the physical features of the studied liquid and 
have been explicitly calculated within MCT for several 
models (see for example for hard spheres [§3|). The a- 
relaxation time scale r is a temperature dependent pa- 
rameter which, according to MCT, scales as the inverse 
of diffusivity 



-co 



IT - T r 



(15) 



The exponent b is fixed by the value of 7 p5| . The 
Kohlrausch- Williams- Watts stretched exponential form 



4>(t)=A K e w * 



-fx 



(16) 



is often used to empirically fit the last part of the <j>(t) 
decay. 

Since all correlation functions obey the time- 
temperature superposition principle, the parameters of 
the two functional forms are calculated by fitting the 
T = 275 K correlation function in the time interval 
30 < t < 600 ps for the Von Schweidler law and in 



the time interval t > 20 ps for the stretched exponential 
form. Similar results are obtained by fitting the T = 266 
K data. The fit of the data according to Eq. [l4| has been 
performed by constraining the exponent b to two differ- 
ent values: the value b = 0.64, obtained theoretically in 



Sec. [II for both approximations and the value b = 0.77 



consistent with the exponent 7 extracted from the analy- 
sis of the diffusivity data. The fitting parameters f c does 
not depend on the value of b. Except for a g-independent 
multiplicative factor, also the g-dependence of hi is in- 
dependent of the choice of b. 

Fig. [5] shows the g-dependence of the fitting parame- 
ters according to Eq. |lj. In the case of the self corre- 
lation function, the non-ergodicity parameter f c (q) pro- 
vides a description, in g-space, of the confining cage. The 
width of the COM cage is larger than the site one, sup- 
porting the view that the site experience an additional 
derealization associated to the hindered librational mo- 
tions. Moreover, we note that the site non-ergodicity pa- 
rameter is not described by a simple gaussian in g-space. 
The g-dependence of hi and /12 follows the oscillation of 
the corresponding static structure factor. 

Fig. |^ shows the fitting parameters according to Eq. 
|l6| . At small wavevectors the extremely long decay of 
the correlation function, controlled by the diffusion of the 
molecules over distances of the order of g _1 does not al- 
low an unbiased determination of the fitting parameters. 
From the wavevector range where the fitting parameters 
are reliable, one can notice that, as commonly found, 
Ak mimics f qi that (3 goes to 1 at small q and that the 
product Dq 2 r has a weak g-dependence. 

We also report the rotational dynamics of the LW 
molecule along the symmetry axis. Fig. [7](top) shows the 
time dependence of the first five Legendre polynomials. 
We also show the corresponding fit with the the stretched 
exponential and the Von Schweidler law. For all angular 
correlators, the time-temperature superposition principle 
holds beautifully, as shown in Fig. 0(bottom). Since the 
quality of the data is very good, even tiny discrepancies 
would be visible. For completeness, we report also the 
fitting parameters according to the functional forms of 
Eqs. 
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in Figs. |8|-|9| The monotonic behavior of the 
fitting parameters in I suggests that for the LW potential, 
the scenario for the rotational dynamics is of the strong 
hindrance type H), as found also for the case of SPC/E 
water ||@. 

The theoretical prediction of the wavevector depen- 
dence of the quantities reported in Figs. |||||| an d % is 
one of the major challenge to the recent proposed molec- 
ular mode coupling theories ^?],|32|,[33| . 



C. Collective Dynamics 

The collective density fluctuations are particularly rel- 
evant in the description of liquids dynamics . The non- 
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monotonic wavevector dependence of the characteristic 
times provides indication of the different timescale of the 
structural modes of the system and their connection with 
the system characteristic lengths. 

We start discussing the site-site correlation functions, 
which can be compared with the corresponding corre- 
lation functions calculated according to MCT, as dis- 
cussed in Sec. HI. In comparing with the Von Schwci- 



dler and the stretched exponential parameterizations, the 
solutions of the MCT equations (with and without the 
C3 contribution), have been treated similarly to the MD 
data. Fig. [1^ shows the q-dependence of the stretched 
exponential parameters, Ak, Pk and tk- We note that 
the theory captures the g-dependence of the a-relaxation 
phenomenon close to the maximum of the site-site struc- 
ture factor. The large wavevector region is also described 
in a reasonable way. The theory fails in describing the 
non-ergodicity parameters in the region where the COM 
structure factor has his peak position. The slowest modes 
in the system are indeed located around 9 nm" 1 {S COM 
peak position) and around 15 nm" 1 (S s peak position). 
The role of the center of mass static correlation in con- 
trolling the molecular dynamics has been observed in re- 
cent experimental studies on supercooled molecular liq- 
uids @. 

The slowing down of the dynamics in the very low-g 
region, below the region of influence of the COM, is also 
predicted in a reasonable way by the site MCT developed 
in Sec. 



Ill 



The stretching exponent parameter (5k is the 
one which seems to suffer more by the chosen approxi- 
mations, which neglects the site-site intramolecular con- 
strains. Indeed, the theoretical estimate is always higher 
than the MD results. The prediction of a higher value 
of (3 is consistent with the reduced number of vertices 
included in the chosen MCT approach. A more elabo- 
rate approach would require as input a large number of 
angular static structure factor correlators |49[1 (which for 
the present molecule would make the theory difficult to 
solve numerically) or the site-site partial static structure 
factors, discriminating at least between the central and 
the two external sites. The development of the molecu- 
lar and site-site MCT approaches is taking place in the 
present days and we look forward applying them to the 
LW model. 

To substantiate the hypothesis that the disagreement 
between theory and MD data in the region around 9 
nmT 1 arises from the approximation used, which ne- 
glects the constrains among sites due to the intramolec- 
ular bonding, we compare in Fig. ^ the de Gennes 
predictions for the COM and site correlation times with 
the corresponding MD data (from Fig. |T^ ). It is clearly 
seen that the dynamics in the wavevector region around 
9 nm -1 is indeed controlled by the COM static structure 
factor @. 

We next compare in Fig. [l(] and |ll| the prediction of 
MCT for the non ergodicity parameter and for the am- 



plitude of the Von Schweidler law with the MD results. 
Again, we notice that the theoretical predictions satisfac- 
torily describe the MD data except for the region where 
the COM structure factor has its maximum. Including 
the C3 contribution improves the agreement in the hydro- 
dynamic limit. Around the peak position of S s (q), the 
agreement is very good, considering that neither fitting 
nor scaling parameters are involved in the comparison. In 
this wavevector region, the long-time dependence of the 
collective site-site correlation function F s (q,t) is rather 
well predicted by the theory, as shown in Fig. OL 

We now turn to the COM correlators. These functions 
cannot be compared with the theoretical predictions re- 
ported in Sec. Ill, since this approximation neglects the 
molecular constrains. But the g-dependence of the COM 
correlators is very relevant for comparisons with future 
molecular MCT predictions. Figs. 15| and 16 show the 
parameterization of the COM collective a-relaxation us- 
ing Eqs. |1] and |l4|. The rather featureless shape of 
the COM structure, which as seen in Fig. |l| shows only 
a broad peak around 9 nm -1 , carries on to the the q- 
dependence of all fitting parameters. 



V. CONCLUSIONS 

The present study of the simple molecular model, in- 
troduced by Lewis and Wahnstrom to mimic the behav- 
ior of ort/10-tcrphcnyl, provides a clean and simple refer- 
ence model for theoretical descriptions of the dynamics 
of molecular liquids in supercooled states. 

The simple MCT analysis reported in this article and 
compared in detail with the MD data shows that already 
the MCT equations for atomic liquids are able to provide 
a rather accurate description of the site dynamics, except 
for the region around 9 nmT 1 where the COM dynamics 
is dominant. We have provided evidence that, in this re- 
gion, the center of mass fully controls the slow-dynamics 
in the liquid Taking into account the triplet corre- 
lation function C3 is shown to improve significantly only 
the small wavevector region. This is not surprising, since 
the usually employed approximation C3 = fails badly 
at small wavevectors |5l}| . While this is not a major is- 
sue in the intermediate wavevector range (i.e. for the 
caging physics) it will become an important element in 
the attempt to extend (toward long wavelengths) the the- 
oretical predictions. 

The present theoretical work, which can be considered 
a zeroth order approximation to the molecular descrip- 
tion, should be followed by an accurate site-site MCT 
approach |?2| or by a full molecular MCT approach |}4| . 
Theoretical development along these two lines are taking 
place at a fast pace and in the near future accurate the- 
oretical predictions should be available for the quantities 
studied in this article. For the present time, we can only 
state that the characteristic MCT scenario appears able 
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to describe the dynamics of this simple, yet complete, 
molecular model. 
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FIG. 1. S s (q) and S, 
data for T = 255K, T 



com(i)- F° r clarity reasons, only 
= 281^ and T = 346K are shown. 



Note that S COM (q) is temperature independent, while for the 
S s (q), the height of the first peak increases on cooling. The 
inset shows the T-dependence of the structure factor at q — 0. 
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FIG. 2. Temperature dependence of the diffusion coeffi- 
cient, evaluated from the long time limit of the mean square 
displacement. The inset show the same data as a function of 
1/T. 
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3. Site self-correlation function for all studied tem- 



peratures at q = 15 nm (the location of the first maximum 
of the site structure factor). 
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FIG. 4. a-scaling for the COM self-correlation function for 
different T at q = 13 nm -1 — the position of the first min- 
imum in ScoM(q)- The main plot shows, from rigth to left 
T = 346, 318, 305, 291, 281, 275K. The inset shows the region 
around the plateau for T = 275, 266, 260, 255^. 
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FIG. 5. Wavevector dependence of the fitting parameters 
according to the Von Schweidler law (Eq. |l4| ) for the self cor- 
relation functions F° elf (q,t) and F** f M (q,t). T = 275K. 
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FIG. 6. Wavevector dependence of the fitting parame- 
ters according to the stretches exponential form (Eq. 
iorFr lf (q, t) and FfJ M (q, t).T = 275X. 
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FIG. 7. (top) Time dependence of the first five Legendre 
polynomials at T = 275K. (bottom) Time- Temperature su- 
perposition for the second-order Legendre polynomium. 
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FIG. 8. Fitting parameters according to the Von Schwei- 
dler law of the first five Legendre polynomials 
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FIG. 10. Wavevector dependence of the non-ergodicity fac- 
tor. The MD data has been calculated as fitting parameters 
according to the Von Schweidler law ( |l4| ) The full and dashed 
lines are the prediction of the MCT using the site structure 
factor as input and including or excluding the C3 contribution. 
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FIG. 9. Legendre's order dependence of the fitting param- 
eters according to the stretches exponential form (Eq. [l6]) 
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FIG. 11. Wavevector dependence of the amplitude hi . The 
MD data (symbols) are the result of the fit according to the 
Von Schweidler law (Eq. ^). Lines are MCT predictions. 
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FIG. 12. Site-site a-relaxation dynamics at the peak of the 
structure factor at T = 275K. Symbols are MD results, lines 
are MCT predictions. 
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FIG. 13. Wavevector dependence of the fitting parameters 
according to the stretches exponential form (Eq. |l6| ) for the 
normalized F a (q,t). Symbols are MD results. The full and 
dashed lines are the prediction of the MCT using the site 
structure factor as input and including or excluding the cz 
contribution. 
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FIG. 14. Comparison of the wavevector dependence of the 
a-relaxation time with the de Gennes approximation for both 
site and COM. 
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FIG. 15. Wavevector dependence of the fitting parameters 
according to the stretches exponential form (Eq. Dm for the 



normalized F„ 



,(q,t). 
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FIG. 16. Wavevector dependence of the fitting parameters 
according to Eq. [m] for the normalized F COM (q, t). 

TABLE I. For each temperature we list the diffusion 
coefficient, the density, the averages of the internal en- 
ergy and of the external pressure and the lenght of the 
simulation time after the equilibration. 
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